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The  mechanical  response  of  living  tissue  is  important  to  understanding  the  injury-risk 
associated  with  impact  events.  Often,  ballistic  gelatin  or  synthetic  materials  are 
developed  to  serve  as  tissue  surrogates  in  mechanical  testing.  Unfortunately,  current 
materials  are  not  optimal  and  present  several  experimental  challenges.  Bulk 
measurement  techniques,  such  as  compression  and  shear  testing  geometries,  do  not 
fully  represent  the  stress  states  and  rate  of  loading  experienced  in  an  actual  impact 
event.  Indentation  testing  induces  deviatoric  stress  states  as  well  as  strain  rates  not 
typically  available  to  bulk  measurement  equipment.  In  this  work,  a  ballistic  gelatin  and 
two  styrene-isoprene  triblock  copolymer  gels  are  tested  and  compared  using  both 
macroscale  and  microscale  measurements.  A  methodology  is  presented  to  conduct 
instrumented  indentation  experiments  on  materials  with  a  modulus  far  below  1  MPa. 
The  synthetic  triblock  copolymer  gels  were  much  easier  to  test  than  the  ballistic 
gelatin.  Compared  to  ballistic  gelatin,  both  copolymer  gels  were  found  to  have  a 
greater  degree  of  thermal  stability.  All  of  the  materials  exhibit  strain-rate  dependence, 
although  the  magnitude  of  dependence  was  a  function  of  the  loading  rate  and  testing 
method. 


I.  INTRODUCTION 

Blunt  trauma  may  have  significant  detrimental  effects 
on  thoracic  organs,  and  understanding  these  effects  is 
critical  for  the  design  of  protective  equipment  used  in 
automotive,  body  armor,  and  sports  applications,  among 
others.  Improved  understanding  of  energy  dissipation 
and  damage  assessment  is  needed  to  model  interactions 
between  human  soft  tissue,  personnel  protective  equip¬ 
ment,  and  the  impact  source.  Experiments  using  cadav¬ 
ers,  instrumented  dummies,  and  synthetic  biomimetic 
materials  measure  and  model  effects  of  extreme  pressure 
and  velocity  on  soft  tissue.  As  a  result,  standard  tests  and 
surrogate  materials  have  been  developed  to  evaluate  and 
rank  protective  equipment.1 
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Recent  work  to  better  understand  how  impact  forces 
are  transmitted  through  soft  tissue  have  led  to  the  fabri¬ 
cation  of  instrumented  surrogate  torsos  composed  of 
modified  ballistic  gelatin  (BG)2  or  modified  silicones.3 
To  an  extent,  these  testing  platforms  mimic  the  mechani¬ 
cal  response  of  soft  tissue  and  are  often  instrumented  to 
gather  pressure  and  velocity  data  during  impact.  This 
data  is  further  utilized  to  validate  finite  element  models, 
which  help  to  describe  the  response  of  the  soft  tissue. 
Often  lacking  in  these  efforts  is  the  proper  match  of  the 
surrogate  mechanical  properties  to  soft  tissue  mechanical 
properties.  A  closer  pairing  of  the  acoustic  and  mechani¬ 
cal  properties  of  a  well-defined  synthetic  material  to 
natural  tissue  can  improve  the  data  collected  from  the 
surrogates,  ultimately  improving  the  ability  to  model  im¬ 
pact  damage. 

Naturally  derived  BG  has  been  most  widely  used  as  a 
tissue  surrogate  to  develop  a  wound  profile  that  provides 
limited  injury  analysis.4,5  For  over  40  years,  BG  has 
been  a  standard  test  medium  for  evaluating  the  effects 
of  ballistics  and  firearms  on  soft  tissue.6,7  BG  is  a 
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thermoreversible  material  formed  from  polypeptide 
chains  dissolved  in  water.  It  is  crosslinked  by  collagen¬ 
like  strands  that  aggregate  to  form  a  triple  helix  structure. 
Previous  research  has  shown  that  the  elastic  modulus  of 
soft  tissue  varies  from  25  to  300  kPa8  with  the  modu¬ 
lus  of  BG  being  approximately  100-150  kPa.9,10  It  is 
readily  available,  cheap,  and  easy  to  process.  However, 
challenges  in  using  BG  as  a  standard  test  medium  are 
significant.  The  mechanical  properties,  for  example,  are 
related  to  the  concentration  of  water  in  the  gelatin,  and  so 
water  evaporation  dramatically  alters  the  mechanical 
properties.  Ideally,  model  gels  would  be  tested  at  37  °C, 
normal  body  temperature,  but  gelatin  provides  usable  tis¬ 
sue  mimetic  properties  only  within  a  narrow  temperature 
range  of  3-30  °C.  Naturally  derived  gelatin  is  polydis- 
perse  and  can  have  a  variable  molecular-weight  distri¬ 
bution  that  can  lead  to  inconsistent  gelatin  properties 
from  batch  to  batch.  Finally,  while  a  properly  formed 
natural  gelatin  is  a  homogenous  material,  procedural  dif¬ 
ferences  exist  with  regard  to  gelatin  processing  and  test¬ 
ing.2,7-10  These  procedural  differences  lead  to  gelatins  of 
differing  mechanical  properties  and  make  comparisons 
between  manufacturers,  testing  labs,  and  researchers  dif¬ 
ficult. 

Any  surrogate  material  selected  to  mitigate  these  is¬ 
sues  and  replace  BG  should  possess  a  well-defined  mo¬ 
lecular  weight  or  crosslink  density  to  ensure  identical 
gels  from  batch  to  batch.  The  material  should  exhibit 
stable  mechanical  properties  as  a  function  of  time  and 
temperature  relevant  to  testing  conditions.  The  material 
should  also  mimic  natural  tissue  mechanical  properties  as 
closely  as  possible.  Ultimately,  this  is  a  challenge  be¬ 
cause  biological  tissue  is  a  non-isotropic  material.  Tri¬ 
block  copolymers  possessing  an  A-B-A  polymer  struc¬ 
ture,  dissolved  in  a  solvent  selective  for  block  B ,  are 
potential  candidates  as  surrogate  materials.  The  elastic 
modulus  of  these  gels  is  primarily  determined  by  the 
molecular  weight  of  the  B  block  (molecular  weight  be¬ 
tween  entanglements)  and  the  polymer  concentration.  By 
adjusting  the  concentration  of  polymer  in  the  gel  or  the 
molecular  weight  of  the  B  blocks,  the  elastic  modulus  of 
the  gels  can  be  tailored  to  mimic  the  elastic  modulus  of 
different  types  of  soft  tissue. 

The  time  and  strain  dependent  mechanical  properties 
of  these  gels  are  determined  by  the  temperature- 
dependent  solubility  of  the  A  block  in  the  selective  sol¬ 
vent.  At  elevated  temperatures,  where  both  the  A  and  B 
chains  are  dissolved,  a  freely  flowing  solution  exists.  As 
the  temperature  is  lowered  through  the  order-disorder 
temperature,  Todt,  the  A  block  chains  precipitate  and 
aggregate  giving  rise  to  the  network  structure  responsible 
for  the  gel  elasticity.  In  this  state,  however,  the  A  block 
aggregates  are  still  swollen  with  solvent.  The  network 
structure  is  viscoelastic  because  the  A  block  chains  can 
exchange  between  aggregates  to  minimize  an  applied 


stress.11,12  As  the  temperature  is  lowered  further  from  the 
^ODT’  the  energy  required  to  pull  chains  from  the  A  ag¬ 
gregates  increases  as  the  A  domains  become  less  soluble 
and  less  swollen.  If  the  temperature  is  further  lowered 
below  the  glass  transition  temperature  of  the  A  blocks, 
the  gel  becomes  highly  elastic  and  exhibits  long  relaxa¬ 
tion  times  characteristic  of  the  glassy  nature  of  the  ag¬ 
gregates.  The  mechanical  properties  of  a  variety  of  ABA 
type  gels  have  been  extensively  evaluated  in  the  litera¬ 
ture.13-18  A  triblock  copolymer  gel  serves  as  an  excellent 
starting  point  for  a  surrogate  material  because  these  gels 
exhibit  predictable  mechanical  behavior  over  a  known 
temperature  range,  the  mechanical  properties  are  rela¬ 
tively  stable  over  time,  and  the  gel  mechanical  properties 
may  be  tuned  by  changing  solvent  selectivity,  triblock 
copolymer  composition,  and  temperature. 

Techniques  such  as  quasi-static  and  dynamic  compres¬ 
sion  and  torsion  tests  can  be  used  to  evaluate  the  me¬ 
chanical  properties  of  gels  on  the  macro-scale,  and  in¬ 
strumented  indentation19  can  be  utilized  to  measure  the 
mechanical  properties  of  these  materials  at  the  micro¬ 
scale.  The  size  scale  of  the  indenter  is  adequate  for  in¬ 
vestigating  heterogeneity  on  the  order  of  microns.20,21 
Instrumented  indentation  has  been  used  to  characterize 
quasi-static  and  dynamic  properties  of  thermoset  poly¬ 
mers  and  elastomers.22,23  There  have  been  numerous  ef¬ 
forts  to  model  the  viscoelastic  response  of  polymers  to 
indentation  based  on  the  use  of  Kelvin  and  Voigt  ele¬ 
ments  to  develop  constitutive  relations.24-27  In  the  case 
of  polymeric  materials,  the  tip-sample  interaction  is  of 
critical  importance28  because  polymers  exhibit  strain  rate 
sensitivity  and  time-dependent  mechanical  phenomena 
such  as  creep  and  stress  relaxation. 

The  use  of  instrumented  indentation  for  evaluating  tis¬ 
sue  surrogates  is  advantageous  because  the  surface  com¬ 
pression  loading  is  similar  to  that  of  a  blunt  impact  into 
soft  tissue.  Also,  current  instrumented  indentation  testing 
machines  permit  the  exploration  of  strain  rates,  stress 
states,  and  resolution  of  measurement  not  achievable  uti¬ 
lizing  conventional  mechanical  tests,  although  the  cur¬ 
rent  methodologies  for  modeling  tip-sample  interactions 
are  limited.29,30  The  theoretical  capabilities  of  these  in¬ 
dentation  systems  include  the  resolution  of  displace¬ 
ments  below  1  nm  and  forces  below  1  p,N,  with  meas¬ 
urements  on  materials  possessing  an  elastic  modulus 
above  1  MPa  readily  achieved.  Significant  challenges 
exist  when  conducting  measurements  on  materials  with 
an  elastic  modulus  below  1  MPa  using  contact  diameters 
of  just  a  few  pan.  Adhesion  between  the  indenter  tip  and 
sample  surface  affect  the  governing  displacement-based 
Hertzian  contact  area  approximation  used  to  determine 
the  contact  area  that  in  turn  affects  the  mechanical  prop¬ 
erty  calculations. 

In  this  work,  the  mechanical  behavior  of  a  traditional 
BG  is  compared  to  two  different  physically  associating 
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gels.  A  methodology  is  described  to  conduct  surface  me¬ 
chanical  property  measurements,  using  instrumented  in¬ 
dentation,  on  soft  gel  systems  (E  <  1  MPa).  Second,  a 
comparison  is  made  using  both  bulk  and  surface  meas¬ 
urement  techniques  to  highlight  the  similarities  and  dif¬ 
ferences  between  the  materials  at  differing  length  scales 
and  frequency  regimes. 

II.  THEORY 

The  governing  equations  for  determining  the  modulus 
of  elasticity  from  compression  and  torsion  tests  are 
readily  available  in  the  literature.31  Thus,  the  discussion 
that  proceeds  is  focused  on  the  contact  mechanics  re¬ 
quired  to  measure  the  mechanical  properties  of  gels 
through  instrumented  indentation  utilizing  a  cylindrical 
flat  punch  in  an  axisymmetric  geometry.  The  cylindrical 
punch  was  chosen  because  the  uncertainty  in  the  contact 
area  is  not  significant  compared  to  other  tip  geometries, 
leading  to  smaller  error  in  measured  mechanical  proper¬ 
ties. 

When  analyzing  bulk  compression  test  data,  load  and 
displacement  measurements  must  be  converted  into 
stress  and  strain  to  obtain  the  modulus  of  elasticity.  For 
a  cylindrical  specimen,  the  average  engineering  stress  is 
computed  as  the  load  divided  by  the  cross-sectional  con¬ 
tact  area  and  the  engineering  strain  as  the  ratio  of  change 
in  specimen  length  to  original  length.  The  modulus  val¬ 
ues  reported  in  this  work  are  the  ratio  of  the  average 
engineering  stress  to  engineering  strain  in  the  elastic  re¬ 
gime. 

To  convert  indentation  loading  data  into  a  modulus 
value,  a  relation  from  contact  mechanics20  is  used: 

=  L  .  a) 

where  E*  is  the  reduced  modulus  of  the  indented  mate¬ 
rial,  r  is  the  cylindrical  radius  of  the  indenter,  and  h  is  the 
displacement  into  the  surface.  The  reduced  modulus  is 
defined  as 


~l-v2  1  -  v2 

E*  =  l—  +  — 

where  v  and  E  are  the  Poisson’s  ratio  and  elastic  modulus 
of  the  indented  material,  respectively,  and  the  subscript  i 
refers  to  properties  of  the  indenter.  If  the  indenter  tip  is 
taken  to  be  perfectly  rigid,  and  the  material  is  assumed  to 
be  purely  elastic,  isotropic  and  incompressible,  the  in¬ 
dented  material  modulus  becomes 


3  P 
8  rh 


(3) 


For  time-dependent  materials,  the  modulus  measure¬ 
ment  is  represented  by  a  complex  modulus,  or  the  vector 


sum  of  the  storage  and  loss  modulus.  Note  that  this  equa¬ 
tion  is  valid  only  for  a  cylindrical  flat  punch.  Relation¬ 
ships  for  other  geometries,  i.e.,  spherical  or  conical  tips, 
will  have  different  relationships  between  modulus,  load, 
and  displacement,  and  will  also  include  a  continuously 
changing  contact  area,  which  is  often  difficult  to  accu¬ 
rately  predict  for  very  soft  materials  without  optically 
measuring  contact  area. 

Traditionally,  storage  and  loss  moduli  are  calculated 
by  considering  the  phase  difference  8  between  an  applied 
oscillating  shear  stress  or  strain,  and  the  delay  in  material 
response  as  a  function  of  time.  In  shear,  the  storage  and 
loss  modulus  are 


G'  =G  cosS 


and 


G"  =  G  sin8 


(4) 

(5) 


respectively,  where  G  is  the  ratio  of  the  applied  shear 
stress  amplitude  to  shear  strain  amplitude.  To  equate 
these  with  uniaxial  storage  and  loss  moduli  values  for  an 
incompressible  elastic  solid,  the  relation 

E  =  3G  ,  (6) 


is  used. 

For  indentation  with  a  flat  punch,  storage  and  loss 
moduli  are  calculated  in  the  following  way.  Assuming 
isotropic,  elastic  and  incompressible  contact,  storage 
modulus  is  a  function  of  the  dynamically  measured  stiff¬ 
ness  S  (or  an  instantaneous  ratio  of  load  to  displacement) 
and  is  given  by 

3S 


This  is  essentially  identical  to  Eq.  (3),  except  that  stiff¬ 
ness  is  measured  dynamically  instead  of  from  the  loading 
or  unloading  curve.  The  loss  modulus  is  a  function  of  the 
measured  damping,  C,  associated  with  the  contact,  along 
with  the  excitation  frequency  a)  and  is  given  by 

3o)C 


These  relations  have  been  shown  to  be  valid  for  dif¬ 
ferent  polymeric  systems.23,32,33 

Mechanical  properties  provide  an  evaluation  of  the 
structural  integrity  of  a  material  as  a  function  of  experi¬ 
mental  variables.  For  a  blunt  trauma  test,  it  is  desirable  to 
develop  a  measure  of  a  material’s  ability  to  dissipate 
impact  energy  as  a  function  of  both  force  and  displace¬ 
ment  rate.  The  energy  per  unit  volume  may  be  used  to 
describe  the  energy  expended  to  deform  the  surface  and 
sub-surface  material.  Energy  per  unit  volume  is  calcu¬ 
lated  from  a  compression  test  using  the  total  area  under¬ 
neath  the  stress-strain  curve  up  to  a  given  strain  value. 
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This  area  will  have  standard  units  of  J/m3  and  represents 
the  total  amount  of  mechanical  energy  absorbed  by  a 
normalized  volume  of  material.  When  considering  in¬ 
dentation  with  a  flat  punch,  the  area  underneath  the 
load-displacement  curve  is  the  energy  only,  with  units  of 
N*  m  or  J.  To  normalize  this  amount  by  a  unit  volume, 
the  total  displaced  volume  is  determined  by 

V  =  it  r2h  .  (9) 

The  relations  in  this  section  are  used  extensively  to 
compute  reported  values  in  this  work. 

III.  EXPERIMENTAL 

Ballistic  gelatin  (BG)  samples  were  made  from  a  20% 
mass  fraction  of  250  bloom  type  A  ordnance  gelatin  dis¬ 
solved  into  40  °C  ultrapure  filtered  water.  The  aqueous 
solution  was  centrifuged  for  15  s  to  remove  solids.  Right 
cylinder  compression  specimens  were  fabricated  by 
pouring  the  solution  into  an  open-faced  aluminum 
mold  and  cooling  to  ambient  conditions.  The  diameter 
and  length  of  the  specimen  geometry  are  50.8  mm  and 
12.7  mm,  respectively,  thus  yielding  a  length  to  diameter 
ratio  of  4. 

Physically  associating  gels  were  made  from  two  tri¬ 
block  copolymers  (Kraton  Polymers,  Houston,  TX)  con¬ 
sisting  of  polystyrene  (PS)  and  polyisoprene  (PI).  Both 
triblock  copolymers  consisted  of  80%  mass  fraction 
triblock  (PS-PI-PS)  and  20%  mass  fraction  diblock 
(PS-PI)  composition.  One  gel  contained  15%  by  mass  PS 
(PS- 15),  and  the  other  gel  contained  30%  by  mass  PS 
(PS-30),  both  used  as  received.  The  polymer  was  mixed 
with  white,  light  mineral  oil  (Mallinckrodt  Chemicals, 
St.  Louis,  MO)  at  a  ratio  of  20%  by  volume  polymer  and 
80%  by  volume  mineral  oil.  The  solution  was  placed  in 
a  nitrogen-purged  vacuum  oven  at  -150  °C  and  fully 
dissolved  over  a  period  of  about  6  h,  being  stirred  every 
hour.  The  melt  was  then  poured  onto  a  flat  surface  where 
it  was  allowed  to  cool  and  gelate,  with  the  PI  chains 
bridging  aggregate  PS  crosslinks.  Both  materials  were 
visually  transparent  and  contained  no  detectable  bubbles 
in  their  final  state. 

The  bulk  mechanical  properties  of  the  BG,  PS-15,  and 
PS-30  were  obtained  using  an  Advanced  Rheometer 
2000  (TA  Instruments,  New  Castle,  DE)  and  a  servo- 
hydraulic  Instron  1331  mechanical  testing  frame  (MTS 
Systems,  Eden  Prairie,  MN).  The  rheometer  has  a  force 
resolution  of  about  10  |xN  and  a  displacement  resolution 
of  about  1  pan,  capable  of  both  compression  and  tension 
testing.  A  4.44  kN  capacity  load  cell  was  used  with  the 
Instron  machine  to  reduce  the  signal  to  noise  ratio  when 
testing  soft  materials.  The  force  resolution  of  the  load 
cell  is  about  0.5  mN,  and  displacement  resolution  is 
about  3  pm. 

BG  was  tested  at  10  ±  0.1  °C  in  compression  using  the 
Instron  machine  at  controlled  strain  rates  of  0.001  and 


0.01  s_1.  Twenty  different  samples  comprised  the  data 
set.  Compression  tests  for  PS- 15  and  PS-30  were  per¬ 
formed  on  the  rheometer  at  room  temperature  using  12- 
mm-diameter  cylindrical  specimens  that  had  a  thickness 
of  approximately  3  mm.  Tests  were  displacement  con¬ 
trolled  at  rates  of  5,  10,  50,  and  100  pm/s,  until  a  total  of 
500  pm  displacement  occurred,  yielding  equivalent 
strain  rates  of  approximately  0.0017,  0.0033,  0.0167,  and 
0.0333  s“\  respectively.  In  compression  experiments  for 
BG,  PS- 15,  and  PS-30,  the  loading  surfaces  of  the  test 
samples  were  well  lubricated  with  mineral  oil  to  ensure  a 
uniaxial  stress  state  and  to  reduce  frictional  effects.  The 
diameter  of  the  loading  plates  was  larger  than  the  speci¬ 
men  diameter  and  was  polished  to  a  mirror  finish  to 
reduce  the  friction  between  the  specimen  and  loading 
platens.  Reported  modulus  values  were  computed  for 
each  compression  test  from  a  least-squares  fit  of  the  ex¬ 
perimental  data.  To  further  investigate  the  mechanical 
properties  of  these  three  materials,  dynamic  tests  were 
performed  on  the  rheometer  within  the  temperature  range 
of  0-50  °C  in  torsion  using  oscillation  frequencies  of  0.1 
and  10  Hz. 

PS- 15  and  PS-30  mechanical  properties  were  meas¬ 
ured  at  the  micro-scale  using  a  TriboScope  depth-sensing 
indenter  (Hysitron,  Inc.,  Minneapolis,  MN).  Load  and 
displacement  are  continuously  measured  during  material 
loading  and  unloading  at  resolutions  in  the  sub-|mN  and 
sub-nm  range,  respectively.  For  all  tests  performed  in 
this  work,  a  250- pan  radius  flat  punch  indenter  was  used. 
The  material  surface  was  located  by  manually  lowering 
the  tip  until  an  adhesive  tensile  force  was  sensed 
(a  negative  force).  The  indenter  was  lowered  until  the 
load  reading  became  zero  (Fig.  1).  Quasistatic  indenta¬ 
tion  tests  to  a  maximum  of  600  nm  were  carried  out  at  1, 
10,  25,  50,  100,  150,  and  200  |mN/s  controlled  loading 
rates.  Unloading  rates  were  approximately  the  same  for 
each  test,  and  holding  times  were  less  than  four  seconds. 
Data  points  were  collected  using  active  load  control, 
meaning  the  spring  force  was  subtracted  from  the  meas¬ 
ured  load  in  real  time  to  achieve  true  loading  rates.  Dy¬ 
namic  tests  were  performed  at  depths  of  approximately 
500  nm  into  the  sample  surface  with  oscillation  ampli¬ 
tudes  of  about  15  nm.  Oscillation  frequency  varied  from 
10  to  200  Hz.  Indentation  data  are  not  reported  on  BG 
due  to  non-repeatability  of  tests.  It  was  experimentally 
difficult  to  obtain  repeatable  tests  since  the  properties  of 
BG  change  so  drastically  with  temperature  (see  Fig.  2) 
and  relative  humidity.  The  mechanical  and  geometrical 
stabilities  both  were  affected  by  inhomogeneity  of  the 
gelatin  and  evaporation  of  water  from  the  surface. 

IV.  RESULTS  AND  DISCUSSION 

Bulk  compression  measurements  at  -10  °C  indicate 
BG  modulus  values  of  96  ±  12  kPa  at  a  strain  rate  of 
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FIG. 


tension  on 
indenter 

1.  Evolution  of  measured  forces  between  the  indenter  and  gel  in  determining  the  initial  point  of  contact  and  during  loading. 


0.001  s-1  and  124  ±  4  kPa  at  a  strain  rate  of  0.01  s-1  to 
ultimate  strains  of  about  1.  The  dependence  of  modulus 
on  strain  rate  suggests  that  BG  has  a  significant  vis¬ 
coelastic  component  to  its  mechanical  behavior.  Dy¬ 
namic  measurements  from  rheometer  experiments  con¬ 
firm  this  as  well.  At  0.1  and  10  Hz  excitation  frequencies, 
the  corresponding  loss  modulus  values  were  calculated  to 
be  approximately  20%  of  the  storage  modulus  at  10  °C 
(Fig.  2).  Storage  modulus  was  about  150  kPa  while  the 
loss  modulus  was  about  30  kPa.  In  theory,  the  complex 
modulus  (or  just  “modulus”  as  it  is  referred  to  here)  is  the 
vector  sum  of  the  storage  and  loss  modulus,  which  are 
perpendicular  to  each  other,  or  about  153  kPa  in  this  case. 
The  modulus  measured  by  compression  should  be,  and 
is,  equal  to  or  less  than  this  value.  A  major  weakness  of 
BG,  as  can  be  seen  from  Fig.  2,  is  the  loss  of  the  gelatin 
structure  between  25  and  30  °C  over  a  period  of  time, 
making  testing  at  room  temperature  unreliable. 


FIG.  2.  Storage  and  loss  modulus  properties  for  BG  as  measured 
dynamically  at  0.1  and  10  Hz  with  the  rheometer  from  0  to  30  °C. 
Properties  vary  greatly  as  a  function  of  temperature 


Gel  compression  tests  at  room  temperature  for  both 
PS- 15  and  PS-30  showed  virtually  no  rate  dependence  in 
the  considered  loading  range  of  5-100  pim/s.  Figure  3 
illustrates  the  averaged  modulus  data  along  with  standard 
deviations  for  ten  different  runs.  For  all  tests,  the  average 
modulus  value  for  PS- 15  is  lower  than  that  for  PS-30; 
however,  for  two  of  the  loading  rate  conditions,  the  val¬ 
ues  are  within  the  standard  deviation  of  the  measure¬ 
ments.  As  compared  to  BG,  the  synthetic  gel  modulus 
values  are  significantly  lower  (by  a  factor  of  4-5).  Dy¬ 
namic  measurements,  shown  in  Fig.  4,  agree  with  the 
compression  test  data  trends  at  room  temperature.  The 
complex  modulus  values  of  PS- 15  and  PS-30  are  ap¬ 
proximately  18  and  30  kPa,  respectively.  It  is  unclear 
why  such  a  detectable  modulus  difference  exists  using 
the  dynamic  method  as  compared  to  quasistatic  uni¬ 
axial  compression,  although  the  two  techniques  agree 
better  for  the  PS  gels  than  for  the  BG.  This  could  be  due 
to  the  fact  that  the  compression  samples  had  lubricated 


FIG.  3.  Compression  test  data  for  PS-15  and  PS-30  at  loading  rates  of 
5,  10,  50,  and  100  |mm/s.  Error  bars  show  one  standard  deviation. 
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FIG.  4.  Storage  and  loss  modulus  properties  for  (a)  PS- 15  and  (b) 
PS-30  measured  dynamically  at  0.1  and  10  Hz  with  the  rheometer  in 
the  temperature  range  of  0-50  °C. 


boundaries  while  the  dynamically  measured  samples  did 
not.  Dynamic  measurement  may  also  cause  a  localized 
heating  effect  in  the  sample,  inducing  a  lack  of  thermal 
equilibrium  in  the  material.  Unlike  BG,  PS- 15,  and  PS- 
30  have  fairly  stable  mechanical  properties  in  the  range 
of  0-50  °C,  which  make  their  properties  easier  to  assess 
at  room  temperature  compared  to  BG. 

Indentation  measurements  were  conducted  to  compare 
microscale  mechanical  properties  to  those  measured  at 
the  macroscale.  Quasistatic  indentation  tests  on  PS- 15 
and  PS-30  yielded  a  rate-dependent  hardening  response 
that  was  not  seen  with  the  compression  tests.  For  both 
gels,  a  higher  loading  rate  yielded  a  higher  load  at  the 
same  displacement  into  the  surface.  An  example  of  this 
effect  is  seen  in  Fig.  5.  Only  the  linear  portions  of  the 
loading  curves  in  Fig.  5  were  considered  for  modulus 
measurements  (typically  from  approximately  100-600  nm 


Displacement  Into  Surface  (nm) 

FIG.  5.  Loading  and  unloading  curves  for  loading  rates  of  1,  10,  25, 
50,  100,  150,  and  200  |ulN/s  on  PS-15,  illustrating  the  time-dependent 
response  of  the  material.  Results  were  similar  for  PS-30. 

indentation  depth).  The  nonlinear  portions  of  the  loading 
curves  are  believed  to  be  due  to  the  inertia  of  the  indenter 
fixture  and  artifacts  from  the  closed-loop  control  system 
and  thus  are  discarded  from  the  analysis.  Because  these 
measurements  are  based  on  Eq.  (3),  the  measured  modu¬ 
lus  represents  a  vector  sum  of  the  storage  and  loss  modu¬ 
lus.  If  the  assumptions  made  are  correct,  the  P  versus  h 
plot  should  be  linear  for  a  fiat  punch.  Therefore,  the  slope 
of  the  line  is  critical  to  determining  the  measured  modu¬ 
lus,  while  the  origin  of  the  fit  will  not  influence  the 
measured  modulus.  This  is  not  the  case  with  other  in¬ 
denter  geometries  such  as  spherical,  Berkovich,  conical, 
or  cube  corner. 

Modulus  values  for  PS- 15  and  PS-30  as  a  function  of 
loading  rate  are  reported  in  Fig.  6.  Interestingly,  the 
modulus  of  PS- 15  is  lower  than  that  for  PS-30  at  loading 
rates  below  50  piN/s,  and  above  this  rate,  the  modulus 
values  are  within  the  standard  deviation  of  each  other  (~1 
kPa).  It  is  possible  that  at  such  high  loading  rates,  this 
behavior  is  an  artifact  of  the  machine’s  collection  capa¬ 
bilities,  or  it  could  point  to  the  fact  that  the  mobility  for 
each  material  becomes  similar  at  higher  loading  rates. 
For  stress  rate  comparison,  the  stress  rate  during  inden¬ 
tation  loading  at  50  piN/s  was  similar  to  the  rheometer 
loading  rate  of  10  pim/s.  Unfortunately,  this  comparison 
is  not  valid  for  strain  rate  because  the  indentation  strain 
rate  (loading  rate  divided  by  load)  continually  changes 
with  indentation  depth  during  a  constant  loading  rate  test. 

The  major  noticeable  difference  between  the  quasi¬ 
static  compression  data  and  the  indentation  data  for  PS- 
15  and  PS-30  is  the  evidence  for  rate  dependence  in  the 
materials.  Although  the  average  stress  in  the  direction  of 
applied  force  during  indentation  is  much  less  than  for 
compression  tests,  the  stress  field  underneath  the  fiat 
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FIG.  6.  Indentation  modulus  of  PS- 15  and  PS-30  as  a  function  of 
loading  rate.  Error  bars  represent  one  standard  deviation  in  the  meas¬ 
urements. 


punch  indenter  is  much  different  from  that  for  uniaxial 
compression.  For  example,  at  about  85%  of  the  radius 
value,  the  stress  underneath  the  tool  is  twice  the  stress 
incurred  at  the  center  of  the  flat  punch.  Uniaxial  com¬ 
pression  assumes  an  evenly  distributed  stress  state  in  the 
direction  of  loading.  However,  at  the  edge  of  the  flat 
punch,  stresses  are  very  large  (theoretically  infinite). 
Therefore,  this  highly  stressed  annular  boundary  area, 
which  is  not  present  in  bulk  compression  tests,  may  con¬ 
tribute  to  the  rate  dependence  of  the  two  materials. 

Likewise,  the  modulus  measured  by  the  indentation 
technique  is  also  higher  than  that  for  compression.  A 
reason  for  the  higher  modulus  may  be  tensile  stress 
buildup  that  contributes  to  frictional  forces,  as  the  incom¬ 
pressible  gel  material  is  stretched  across  the  face  of  the 
indenter.  The  average  storage  modulus  of  both  PS  gels 
was  47  ±  8  kPa,  and  the  loss  modulus  was  2.4  ±1.0  kPa 
when  the  data  were  averaged  over  the  frequency  range  of 
10-200  Hz,  as  measured  by  dynamic  indentation.  These 
averages  yield  complex  modulus  values  only  slightly 
higher  than  those  found  using  the  quasistatic  indentation 
technique,  which  could  be  due  to  heating  effects  at  higher 
frequencies.  However,  the  measured  storage  and  loss 
modulus  varies  quite  a  bit  over  the  frequency  range 
(Fig.  7).  At  10  Hz,  both  materials  had  a  storage  modulus 
of  39  kPa,  which  is  still  higher  than  values  determined 
using  the  rheometer,  but  are  within  close  agreement  of 
the  quasistatic  data.  This  discrepancy  between  indenta¬ 
tion  storage  modulus  and  the  rheometer  values  may  be 
dependent  on  the  calibrated  dynamic  response  of  the  in¬ 
dentation  instrument,  which  has  a  much  greater  stiff¬ 
ness  than  that  of  the  contact.  Therefore,  the  dynamic 
contact  stiffness  and  damping  measurements  may  not  be 
sufficiently  sensitive  to  properly  sense  the  contact. 


Frequency  (Hz) 

FIG.  7.  Storage  and  loss  modulus  of  PS- 15  as  measured  dynamically 
with  the  flat  punch  indenter  over  a  range  of  10-200  Hz.  Results  for 
PS-30  were  similar. 

Similar  to  the  trend  with  the  loss  modulus  measured  via 
the  rheometer,  there  is  a  consistent  increase  in  its  mag¬ 
nitude  as  a  function  of  frequency.  Time-temperature  su¬ 
perposition  measurements  (not  shown)  made  via  the  rhe¬ 
ometer  and  shifted  to  frequencies  above  100  Hz  show  a 
monotonically  increasing  modulus  compared  to  the  de¬ 
creasing  modulus  observed  in  this  range  for  the  dynamic 
indentation  measurements.  As  measured  with  indenta¬ 
tion,  the  loss  modulus  is  0.20  ±  0.03  kPa  at  10  Hz  and 
3.7  ±  0.1  kPa  at  200  Hz  for  both  PS-15  and  PS-30.  The 
loss  modulus  was  0.5%  of  the  storage  modulus  at  10  Hz 
and  16%  at  200  Hz. 

Because  the  materials  in  this  study  are  used  to  examine 
projectile  penetration  characteristics,  it  is  also  useful  to 
consider  how  much  energy  per  unit  volume  is  absorbed 
by  each  material  during  an  impact.  This  will  not  be  a 
perfect  correlation  because  the  magnitude  of  strain  rates 
seen  during  a  ballistic  event  (-1,000-10,000  s-1)  are 
much  higher  than  those  for  indentation  (-0.01-100  s-1), 
but  insight  can  still  be  gained.  In  either  case,  the  higher 
the  impact  energy,  the  greater  the  toughness  of  the  ma¬ 
terial  required  to  absorb  the  projectile  without  critical 
failure.  The  absorbed  energy  will  change  depending  on 
the  displacement  during  the  compression  test,  the  sample 
thickness,  or  displacement  into  the  surface  for  an  inden¬ 
tation  test.  For  BG,  the  energy  per  unit  volume  as  a 
function  of  displacement  during  compression  for  0.001 
and  0.01  s-1  strain  rates  is  compared  to  test  data  on  PS-15 
and  PS-30  [Fig.  8(a)].  For  PS- 15  and  PS-30,  there  was  no 
energy  per  unit  volume  dependence  on  strain  rate.  Under 
uniaxial  loading,  the  energy  per  unit  volume  increase  is 
exponential  with  displacement.  For  any  loading  geom¬ 
etry,  the  energy  per  unit  volume  magnitude  as  a  function 
of  displacement  trends  with  the  measured  modulus  value 
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FIG.  8.  Energy  per  unit  volume  as  a  function  of  displacement  for  (a) 
a  compression  test  and  (b)  a  flat  punch  indentation  test.  Compression 
test  data  show  a  comparison  between  energy  absorbed  under  similar 
loading  conditions  for  BG  at  0.001  and  0.01  s-1  and  PS- 15  and  PS-30 
(extrapolated  to  300  |ixm  penetrations).  Energy  absorbed  by  PS-15  in 
(b)  was  similar  to  that  absorbed  by  PS-30,  although  PS-30  was  slightly 
higher. 

only.  Figure  8(b)  shows  a  comparison  between  PS- 15 
and  PS-30  for  the  flat  punch  loading  geometry.  Be¬ 
cause  BG  could  not  be  tested  at  the  microscale,  no  similar 
comparison  was  conducted.  PS-30  absorbs  more  energy  at 
the  various  loading  rates  than  PS-15  (i.e.,  has  a  higher 
modulus),  especially  at  lower  loading  rates.  This  is  possibly 
a  reflection  of  the  higher  force  required  to  remove  indi¬ 
vidual  PS  chains  from  larger  crosslinks  of  this  material. 

Between  the  two  different  loading  geometries  used  in 
this  study  (or  between  any  two  loading  geometries), 
the  energy  per  unit  volume  comparison  for  the  same 
material  is  most  useful.  Approximately  two  hundred 
times  more  energy  is  absorbed  in  the  indentation  tests  at 
300  pirn  (extrapolated)  for  PS- 15  and  PS-30  than  in  the 
compression  tests.  This  comparison  suggests  the  degree 


to  which  a  500- pirn  diameter  flat  punch  penetrator  may 
cause  failure  to  an  impacted  material  compared  to  a  large 
blunt  force.  In  the  future,  different  penetrator  geometries 
will  be  examined  to  determine  the  degree  of  impact  dam¬ 
age  from  variously  shaped  projectiles. 

V.  CONCLUSIONS 

The  modulus  of  BG  was  successfully  measured  at 
10  °C  using  quasistatic  and  dynamic  techniques  on  the 
macroscale.  Two  physically  associating  gels,  PS- 15  and 
PS-30,  based  on  swollen  triblock  copolymers,  were 
found  to  exhibit  quasistatic  properties  within  the  range  of 
BG.  Although  the  modulus  values  of  PS-15  and  PS-30 
were  less  than  BG,  they  were  still  at  the  lower  end  of  the 
modulus  range  reported  for  soft  tissue  (~25  kPa).  Even 
though  only  two  specific  gels  were  studied  in  this  work, 
the  gel  systems  may  be  tailored  to  have  their  modulus 
increased.  Compared  to  BG,  both  PS-15  and  PS-30  were 
found  to  have  a  greater  degree  of  thermal  stability.  BG 
gels  have  a  shelf  life  of  just  a  few  days  at  ambient  con¬ 
ditions.  A  successful  methodology  was  developed  for 
mechanically  characterizing  materials  with  a  modulus 
below  1  MPa  using  instrumented  indentation,  which  ap¬ 
plies  to  a  number  of  other  soft  material  systems.  Both 
macro-  and  microscale  tests  of  PS- 15  and  PS-30  yielded 
similar  modulus  measurements.  Macroscale  measure¬ 
ments  found  little  rate  dependence  on  either  PS- 15  or 
PS-30,  while  microscale  measurements  revealed  rate  de¬ 
pendent  behavior.  The  reason  for  this  difference  was 
thought  to  be  due  to  the  boundary  conditions  present 
during  flat  punch  indentation  as  compared  to  a  compres¬ 
sion  test.  Interpretation  of  the  dynamic  mechanical  prop¬ 
erties  of  the  two  gels  measured  with  instrumented  inden¬ 
tation  and  comparisons  to  the  loading  rate  effects  were 
difficult.  More  extensive  work  on  these  and  other  similar 
material  systems  is  planned  for  the  future. 
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